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A NOTE ON QUANTIZATION IN THE PRESENCE OF GRAVITATIONAL 

SHOCK WAVES 

MARCELO M. DISCONZI 

Abstract. We study the quantization of a free scalar field when the background metric 
satisfies Einstein's equations and develops gravitational shock waves. 



1. Introduction. 



The general notion of "singularity" , although not always precisely defined and having dif- 
ferent meanings in different situations, plays a crucial role in many physical theories. Broadly 
speaking, one expects that singularities will form in regimes where the system undergoes 
extreme dynamic conditions, e.g., turbulence in Fluid Dynamics or gravitational collapse in 
General Relativity (GR). 

One particular notion of singular behavior frequently encountered is that of a shock wave. 
r~~~ Roughly, it corresponds to a discontinuity in the solutions of the equations of motion --in 

particular, it implies that such solutions exist only in a generalized sense. In the case of Ein- 
^-j- \ stein's theory of gravity, shock wave discontinuities happen in the first derivatives of the metric 

_il ■ (see definition 2.8). It follows that the spacetime does not carry a "smooth geometry," with 

O ■ the curvature tensor (which depends on second derivatives of the metric) being meaningfully 

defined only in a distributional sense. 

Although in a different context, physicists have been dealing with singular spacetimes for 
quite a long time, with Hawking's singularity theorems [HE] and its ramifications (Cosmic 
Censorship, Penrose Inequality, etc.) being at the core of such developments. It is widely 
believed that a successful quantum theory of gravity will resolve many, if not all, of the 
difficulties and puzzles that arise in singular backgrounds. Unfortunately, despite much of the 
progress that has been witnessed in the last few decades, such a theory is not yet available. 
This does not mean, however, that we cannot learn something about quantum effects in curved 
spacetimes, singular ones included. In the range where curvature effects cannot be neglected 
but are still far from the Planck scale, the powerful (if yet difficult) machinery of Quantum 
Field Theory (QFT) in curved backgrounds is available to us. To this day, some of the best 
hints of what the long sought quantum theory of gravity might look like come from the study 
of quantum fields over a non-flat background that satisfies Einstein's equations. In fact, the 
ability of reproducing the black hole temperature - - discovered by Hawking via a careful 
application of QFT on a background that undergoes classical gravitational collapse [H] - - is 
often regarded as the first test for any theory attempting to quantize the gravitational field. 
This, of course, is not different than many other instances in Physics where semi-classical 
formulations are useful in providing insight into what the full quantum theory is, leading in 
this way to fruitful directions of inquiry. 

l 
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Furthermore, one should not have to wait until a full-fledged theory of quantum gravity is 
in place in order to understand interesting physical phenomena, which involve quantum effects 
in a background where the singularities are still amenable to a fully classical treatment. This 
paper is a step in this direction. We shall study the quantization of a free scalar field over a 
spacetime where gravitational shock waves are formed. Einstein's equations will have to be 
defined in a suitable weak sense, and some standard arguments, like the construction of prop- 
agators, adapted to this weaker setting. A similar aspiration, namely, devise a quantization 
scheme that can potentially incorporate singular background data, was explored in [AD]. 

The tools we shall employ are, in a sense, not new. Most of the constructions go back to 
the work of Lichnerowicz on tensor distributions [LI] (see [D] for similar constructions and 
generalizations). These have the convenience of being simultaneously suited to the study 
of shock waves, on one hand [L2], and to the quantization of fields in curved -- although 
smooth — spacetimes, on the other hand [L3]. Not surprisingly, many of the arguments here 
presented consist of carefully checking that the results of [L3] carry onto the framework of 
shock waves. As such arguments can be done in charts, our point of view will be purely local. 
Generalizations to a global setting are possible, provided that further conditions are taken 
into account; we briefly comment on this at the end. 

Although this paper focuses uniquely on mathematical aspects, we stress that the study of 
gravitational shock waves in general, and corresponding scenarios where quantum effects might 
become important in particular, has attracted significant attention in the Physics community 
(see e.g. [DPI, DP2, DP3, HS, 'tH, HU, VV] and references therein), hence the importance 
of laying out its mathematical foundations. 

Notation 1.1. (i) vol(g) denotes the volume form of the metric g, and \g\ its determinant 
(in a local coordinate patch); (ii) if Q C M, MaC' manifold, T>k tP (£l) denotes the space of 
p-tensors of class C k , < k < £, with compact support in Q, i.e., test tensors in Q; we shall 
use the given metric to identify covariant and contra-variant tensors, hence referring simply 
to "p-tensors"; (iii) V' kp (Q) is the space of continuous real- valued forms on V kjP (Q), where 
continuity is understood in the sense of the theory of distributions [S]; (iv) we sometimes 
write T> p (Q) when the differentiability is clear from the context, although for the most part 
it will suffice to deal with elements in T>o }P (Q) and its dual; (v) (v,u) means u G "Z\ p (fi) 
evaluated at v G V kp (Q), whereas (ui,u 2 ) g denotes the inner product between the p-tensors 
u\ and ii2, although, when the distributions arise from locally integrable tensors, we naturally 
identify (■, ■) and (-, -) 9 ; (vi) we shall denote by T" the tensor distribution defined by a locally 
integrable tensor T; (vii)V will denote covariant differentiation with respect to the metric g. 
(viii) [ck | -y | /?] indicates anti-symmetrization in the indices a and /3. 

Hypotheses in the form 'Assumption X.X" (e.g., Assumption 2.1) are assumed throughout, 
therefore we do not state them in the theorems and definitions. 

2. Gravitational shock waves. 

We start recalling some definitions and fixing out notation. Let M be an oriented 4- 
dimensional differentiable manifold of class C 2 and piecewise C 4 , and Q C M a contractible 
open subset. Whenever coordinates are employed, it is implicitly assumed that Q is taken 
small enough as to belong to the domain of a coordinate patch. Let £ be a regular hypersurface 
in Q which is given locally as the level set {tp = 0} of a C 1 function p> : Q — > R; saying that £ 
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is regular then means dip 7^ on E. We suppose that E partitions Q into two disjoint domains 
Q + and Q~ given by {(p > 0} and {tp < 0}, respectively. 

Definition 2.1. Let T be a p-tensor on Q, and assume that T is of class C° over f2 + and Q~. 
T is regularly discontinuous across 1 E if, when </?—>• + (i.e., 9? tends to zero through positive 
values) (resp. tp —¥ 0~), T converges uniformly to ap-tensor T + (resp. T~) defined on E. The 
discontinuity 2 {{T}} of T is the p-tensor on E defined by {{T}} = T + — T~ . 

Let g be Lorentzian metric 3 of class C° and piecewise C 2 on M. We fix Q, tp and E once 
and for all, and henceforth assume for the rest of the paper: 

Assumption 2.1. (i) tp is C 3 on Q + and Q~ , with second and third derivatives regularly 
discontinuous across E; (ii) g is C° on Q, C 2 on Q + and f2~, and has first and second derivatives 
regularly discontinuous across E. 

Unless stated otherwise, it will be assumed in this section that coordinates {x a } 3 =0 are 
adapted to E, what means that x° = tp and <9, is tangent to E, i = 1, 2, 3. From now on, Greek 
indices run from to 3, while Latin indices run from 1 to 3. 

Definition 2.2. In coordinates {x a } 3 Y=0 adapted to E, the metric components g^ will be 
called physical components with respect to E, while go a will be called gauge (or non-physical) 
components with respect to E. A change of coordinates x' a = x' a (x°, . . . , x 3 ) is called a change 
of gravitational gauge with respect to E, or change of gauge for short, if it is the identity on 
E, and the the values of ^ a ^| s as well as of cfo<7y| E are invariant under this change. 

Remark 2.3. Despite the familiar terminology, the reader should not be led to think that E 
is a space-like hypersurface. In fact, for the case of interest in this paper, E will be light-like. 

The necessity of treating g^ and go a differently comes from the well-known fact that GR has 
a gauge freedom due to the action of the diffeomorphism group of M. In fact, in neighborhood 
of E gauge changes take the form 

x' a (x°, x\x 2 , x 3 ) = x a + ^- (d x' a (0, x\x 2 , x 3 ) + r a (x°, x\x 2 , x 3 )) , 

where r a converges, along with its derivatives, uniformly to zero when x° — >• 0. It is therefore 
possible to arrange the terms in parenthesis as to produce or eliminate discontinuities of dog' 0a 
on E — showing that such terms carry no intrinsic physical meaning [L2]. Furthermore, from 
Assumption 2.1 one readily sees that 

{{9i9afi}} = = di{{ gafi }}. (2.1) 

From this and the above discussion, we see that it suffices to focus on d gij in our study of 
discontinuities of the metric across E — see definition 2.8. 

We now turn our attention to the appropriate notion of weak solution for the Einstein's 
equations in the study of gravitational shock waves. For each fixed pair of indices 05, the 
Christoffel symbols T^ s locally define a vector field T Q = Vgg, which in turn, under our 



"regulicrmcnt discontinues a la traversee de" , [L2] . 

We avoid the more familiar notation [•] for the discontinuity of a function in that it may cause confusion 
when written inside an equation. 
Our convention is (+ — ). 
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hypotheses, is locally integrable and hence defines a vector distribution (P^)'. Following 
Lichnerowicz, it is therefore natural to define the curvature tensor distribution as 

7^87,5 — v 7 (r£g 5 )) — v ( 5(r° /37 )) , 

where the covariant derivatives are interpreted in a distributional sense [LI, L2]. To find a 
simple formula for the associated distributional Ricci curvature and relate it to the ordinary 
one, we shall use the following lemma, whose proof is an application of the tools developed in 

[L2]. 

Lemma 2.4. Let T be a p-tensor of class C 1 on Q + and Q~ such that T and d a T are regularly 
discontinuous across E. Then VT is regularly discontinuous across E. Furthermore, T and 
VT define tensor distributions T' and (VT)' such that 

VT'-(VT)' = d<p®5x{{T}}, 

where the covariant derivative VT' of T' is in the sense of distributions, and 5% is the Dirac 
delta on Q with support on E. 

Proof. Let A 1 ^ = T^dx 1 be the connection one form of g on Q. From our hypotheses, we see 
that A^ is C l in Q + and Q~; it is also regularly discontinuous across E. It follows that VT is 
regularly discontinuous across E since d a T and T are so. 

Let x+ (resp. x~) be the function defined a.e. in Q which equals to 1 in Q + (resp. Q~) 
and zero in Q~ (resp. Q + ). Since T and VT are in Lj oc (Q), they define tensor distributions, 
which can be written as 

T' = x+T + X-T, (2.2) 

and 

(VT)' = X +VT + X-VT, (2.3) 

where these equalities are to be understood in the sense of distributions 4 . We can write 
vol(g) = dtp A u. Notice that u depends on ip, but if rj is another 3-form such that vol(g) = 
dip A rj, then t] = u + dtp A a for some 2-form a. Hence, for any test function u 



uu = j uu 
dn+ Jan- 



(where dQ + and dQ are the oriented boundaries so that dQ + = —dQ ) has a well-defined 
value independent of the choice of u. In particular, (5s, u) is given by 



{Sy,,u} = — / uu — uu . 

Jan+ Jan- 

Notice that because u is continuous this gives that <5s is in fact an element of V 00 (Q). Since 
5-£ can be multiplied by locally integrable functions, this same formula also shows that 5% has 
a well-defined action on 1)^^(0,), k = 1,2. 



Equalities among quantities in T>' (fi) are by definition in the distributional sense, so we shall no longer 
write "in the sense of distributions." 



QUANTIZATION 
Clearly V%+ E Z>i(fi), and for any u E £>i(fi) 

(Vx+,«> = (Vx+,m) s = -(x+,divn) 9 = - / ^^ d a (v/-|ffK) vol(g) 

^n+ V~\9\ 

d a (\/ -\g\u a )dx° A dx 1 A dx 2 Adx 3 = - / ^/Htfl" ^ 1 A dx 2 A dx 3 



(p a u a cu = (dip5x,u), 

so that Vx+ = dip 5^- Similarly, one finds Vx~ — —dtpS-^. Therefore, 

V(x + T) = (Vx+) ® T + x+ VT = d^ s ®T + X+ VT. 

Computing a similar expression for V(x + T) and using (2.2) and (2.3) yields the result. □ 

It is seen that the discontinuities of g enter in the expression for the distributional curvature. 
The following lemma is useful to handle them. 

Lemma 2.5. For each of the functions g a p, there exists a scalar distribution d a p such that 

<5e{{<9 7 # Q/ 3}} = (p-yd al 3. 

Proof. Arguing similarly to the proof of Lemma 2.4 we see that 

(Vj5 E , u) = / di(y/-\g\u) dx 1 A dx 2 A dx 3 
Jdn+ 

did (\/— \g\ u) dx° A dx 1 A dx 2 A dx 3 



= 0, 

where we used that u, being a test function, is compactly supported in Q. A similar result 
holds in Q~ and we conclude that Vjfe = (in a more invariant fashion, we can say Vjfe = 
for vector fields X tangent to tp =constant). It follows that 

V£ s = d(pV 5^. 

This combined with (2.1) (and, of course, our assumptions on g) produces the result. □ 

Intuitively, d a p corresponds to the "jumps" that come when we try to differentiate Sx,{{g a /3}} 
across S. It is not difficult to see that the existence of d a p is equivalent to the formulation 
given in [L2], namely, to find, for each g a p, a function b a p defined over £ such that 

{{d 1 g a /3}} = <P y b a p on £ 

and satisfying 5 d a p = S^b a /3- 

Remark 2.6. b a p depends on the choice of equation for S. The functions b a p were initially 
introduced by Lichnerowicz [L2] and are important in determining the conditions of shock for 
gravitational waves. Such conditions will appear below when we construct the propagator for 
the Klein-Gordon equation. 



Notice that, in principle, b a p is denned only on E, but this suffices since 5s is supported on S; alternatively 
b a p can be extended to 51 without affecting the distributional equality d a p = S^b a /3. 
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From Lemma 2.4, Lemma 2.5 and the above considerations it follows that 

where the distribution R' a g 7S , with Ra^s being the curvature tensor, is well-defined because 
of Assumption 2.1, and Hap^s is the tensor on E given by 

—2% a ^s = b a .yipf3(fs - baStppifiy + bps'Pa'Pj ~ b^ip a <p 5 . 

The distributional Ricci tensor is defined as 

n aB = R' aB + 5^n aB , (2.4) 

where 

2H aB = ba^ipp + bp^tpa - b^cpaipp - b aB ip^ip n, (2.5) 

so that lZ a p is formally the trace of the distributional curvature tensor. 

Recall that given a stress-energy tensor T aB (possibly identically zero), one can write Ein- 
stein's equations as 

Ra/3 = K-PafS, (2.6) 

where k is a constant and 

Pa/3 = Ta/3 ~ ~Tg a /3, 

with T the trace of T a p- The last ingredient we need to define the distributional Einstein's 
equations is the regularity of T a p- 

Assumption 2.2. T a p is a given symmetric two-tensor, continuous on Q + and Q + and regu- 
larly discontinuous across E. 

It follows that p a/ 3 shares the same regularity properties of T a p- The reason why we think 
of T a fi as given is that existence of solutions to the distributional Einstein's equations --in 
which case T a p has a determined functional form but depends on the metric and the matter 
fields of the problem — will not be investigated; rather we shall assume we are given a solution 
that defines a background on which fields will be quantized. 

Definition 2.7. The distributional Einstein's equations are defined as 

Definition 2.8. The hypersurface E is called a wave front and is said to define a gravitational 
shock wave if g satisfies (2.7), and the first derivatives of the physical components of g are 
discontinuous across E, and also regularly discontinuous across E. 

It is possible to show that a shock wave E is necessarily light-like 6 ; had it not been light-like, 
and hence characteristic to the reduced Einstein equations, the values of the induced metric 
and its derivatives would uniquely determine a regular solution on both sides of E, preventing 
discontinuities. See [L2] for details. 



Recall that hydrodynamic shock waves have the property of being supersonic before the shock and subsonic 
after the shock, with ordinary waves propagating at the sound speed. Here, the speed of light plays the role 
of the sound speed, hence shocks cannot be "superluminal" before the shock. 
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3. Quantization. 

In this section, we adapt the techniques of [L3] to the weak setting developed above. We 
assume the same hypotheses and notation as before. We shall deal only with a free scalar 
field, although it is very likely such results can be generalized to tensor and spinor fields. 
Not surprisingly, the extension to interacting theories, on the other hand, is expected to pose 
severe difficulties. For the rest of the paper, we let m > be a fixed parameter. For a compact 
set K C Q, we denote by C + (K) (resp. C~(K)) the future (resp. past) of K in the usual 
sense of GR, and the cone of K the set C(K) = C + (K) U C~(K). As our point of view is 
purely local, strictly speaking C + (K), is the future of K within Vt (analogously for C~(K)). 

Proposition 3.1. For each fixed x E Q,, there exist two elementary kernels E^r satisfying 

(n g + m 2 )E± = 5 x , 

where S x is the Dirac delta supported at x. E+ (resp. E~) is unique and has support on C + (x) 
(resp. C~(x)). 

Remark 3.2. In Minkovski space, E^ 1 are the standard advanced and retarded kernels. 

Proof. We give the proof for E + , with the existence of E~ being completely analogous. First, 
notice that the standard formula 

/ uD g vvo\(g) = / ud li (g' JJ/ y/ -\g\d„v) dx° A dx l A dx 2 A dx 3 = / vD g uvo\(g) 
Jo. Jn Jn 

still holds for all u, v E T>2,o(£V), despite g being only continuous on the whole of Q. D g is, 
therefore, formally self-adjoint, hence we seek to find E£ such that for all u E £>2,o(^) 

(E+,(n g + m 2 )u) = (5 x ,u). 

Following Choquet-Bruhat [Bl], we construct a parametrix er+, i.e., a distribution satisfying 7 

(D g + m 2 )at = 5 X - Qt , (3.1) 

where Q+ is an integrable function supported on T x = dC + (x). Indeed, we shall show that 
the fixed point argument employed in [Bl], where it is assumed that the metric is sufficiently 
different iable, still goes through under our assumptions. 

We start by investigating the restrictions that (2.7) imposes on the the quantities 6 Q/ g - 
these are the so-called conditions of shock, originally studied in [L2]. Because g is C 2 over 
Q + and Q~ , and in light of Assumption 2.1, we have that (2.6) is satisfied in Q + U Q~ and, 
moreover, 

K/3 = Kp' aP in tt. 

From this, (2.4), (2.5) and (2.7) we conclude that 

V^IaHVfl = on S. (3.2) 



The terminology "parametrix" was introduced by Hilbert and usually indicates an approximation for the 
fundamental solution of a linear equation. 
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We now claim that <^V M is a well-defined operator along E. In fact, direct computation gives 
that on E 

which vanishes by (3.2). From this, the regularity of tp in Assumption 2.1 and the fact that 
E is light-like [L2], it follows not only that y^V^ is well-defined over E but also that 

<^V^ Q = on E. 

As a consequence, geodesies are well-defined over (and in fact span) E. Hence, even under the 
circumstance F x D E ^ 0, we can set up a system of ordinary differential equations over T x 
that determines a+; this is done exactly as in [Bl], which the reader is refered to for details. 
From equality (3.1), it now follows that 

u{x) = I Qtuvol(g) + (<r+, (D g + m 2 )u). 
Ju 

This can be viewed as an integral equation for u, which can be solved by a Neumann-series 

type of technique as in chapter V of [dR], possibly after shrinking Q. The solution is unique 

and operates linearly and continuously over u, yielding the desired kernel. □ 

As usual, we think of E^r as defining a distribution E ± (x',x) in Q x Q by 

(E ± (x',x),u(x',x))axU = (E^,(x),u(x',x))uxQ- (3.3) 

The lemma below verifies that the standard symmetry properties of these distributions known 
to be true in the smooth setting, continue to hold under our hypothesis. 

Lemma 3.3. Let E ± (x,x') be as above. Then 

E + (x,x') = E~(x',x), 

for x,x' G Q. 

Proof. Fix a test function u and define 

Vi(x') = (E~(x,x'),u(x)). 

Since E~(x') has support on the past of x, V\ is supported on the past of the support of u. 
We have 

(D x > + m 2 )vi(x') = ({D x > + m 2 )E~(x, x'), u(x)) = (5(x, x'), u(x)) = u(x'), 

where □ = D g , \3 X means that derivatives are with respect to the x variable, and we used the 
symmetry of the Dirac delta. 

For sufficiently smooth metrics, Choquet-Bruhat proved [B2] that any solution i; 2 of 

(□ + m 2 )v 2 = u (3.4) 

with support compact towards the future 8 is given by 

v 2 (x') = (E + (x',x),u(x)), (3.5) 



K C f2 is said compact towards the future if C (K)nC + (x) is compact of empty for each x € $7; a similar 
notion applies for compact towards the past. 
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provided that Q is taken sufficiently small. An inspection on her proof shows that the same 
statement still holds in our setting. In fact, the set K = C_(supp(-u)) fl C + (x) is compact (or 
empty). Choose a function z G V 2t0 (Q) equal to 1 in a compact neighborhood of K. If v 2 
satisfies (3.4), we obtain that for any x' G C_(supp(it)), 

(E + (x', x), u(x)) = (E + (x', x), (D x + m 2 )v 2 (x)) = {E+(a/, x), (D x + m 2 )(z(x)v 2 (x))). 

But, on the other hand, if 

v 2 {x') = (E + (x',x),u(x)}, 

then 

(E + (x' , x) , u(x)) = ((D x + m)E + (x', x), z(x)v 2 (x)) = v 2 (x'), 

which shows (3.5). We conclude that v\ = v 2 , from what the result follows. □ 

Following standard convention and terminology, we then define the propagator of D g + m 2 
as the distribution on Q x Q given by 



It is seen that 



In light of Lemma 3.3, 



G(x, x') = E + (x', x) - E-(x', x) 

(D x + m 2 )G(x,x') = 0. 



G{x,x') = -G{x',x). 

Consider the massive wave equation in O, 

(D g + m 2 )u = 0. (3.6) 

Let us denote by D + (K) (resp. D~(K)) the future (resp. past) domain of dependence of an 
achronal set K, and D(K) = D + (K) U D~(K). We next suppose that: 

Assumption 3.1. There exists in Q a kernel Q(x,x'), which is a solution (in x) of (3.6), 
satisfying Q(x,x') = Q(x',x), and such that, for each space-like three surface S C Q and 
x, x' G D(S) fl £1, the following holds: 

G(x,x') = J {Q(x,y)d»Q{x\y) - G(x',y)d»g(x,y)) dAfty), (3.7) 

where dA g is the volume element induced on S by g. 

Remark 3.4. In Minkovski space, Q is the so-called D 1 distribution associated with the 
Pauli- Jordan propagator [BS]. 

As in the usual case of smooth coefficients, for x G D(S) fl Q we have the following formula 
for a solution u of the Cauchy problem of (3.6) [Bl, B2], 

u{x) = j (u(y)d^G(x,y) - G(x,y)d„u(yf) dA^y). 
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This determines the solution u in terms of the Cauchy data q := u| s and p := 9^w| E . The 
pair (q,p) plays the role of a point on phase space. We now define u, depending on u, by 

and introduce the norm 

Define 
and 



G-iQ „_ G + iQ 



2 ' 



, u — iu _ u + iu 
2 ' 2 

Proposition 3.5. The following identity holds 



^^J = {u ,u ) = -{u,u). 



Proof. This follows directly from the above formulas. □ 

We can now carry out the quantization of a free scalar field over the background (£l,g) by 
following the corresponding steps in the quantization of fields defined over a smooth curved 
background, as originally proposed by Lichnerowicz [L3], and extended by other authors, in 
particular Wald [W] (see also [HW]). In order to obtain a well-defined quantization procedure, 
one has to be specific about the structure of the operators involved, indicating their domains, 
self-adjointness properties etc. However, with the propagator G(x,x'), Proposition 3.5, and 
formula (3.7) at hand, these constructions are the same as in [L3] (see also [W]). Hence, for 
the sake of simplicity, we shall restrict ourselves to indicating what the main features of the 
corresponding quantum theory are. 

We have: 

I. The association u^?u defines a linear automorphism J on the space of solutions 9 of (3.6). 
u ± are eigenfunctions of J corresponding to the eigenvalues ±i. The projections u \— > u^ 
correspond exactly to the decomposition of the field into positive and negative energies. 



9 A 2-form u> can also be introduced in the space of solutions, and shown to define a symplctic structure, u) 
and J are compatible and define an (infinite dimensional) Kahler structure. These intrinsic quantities can be 
used to derive an appropriate notion of inner product to carry out the quantization procedure. See [L3] and 
[W] for details. 
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II. Assume from now on that u is operator- valued. Our definitions and formula (3.7) imply 

[m + ,-u + ] = = [u~,u~], 

[u + (x), u"(x')] = —iG + (x,x')id, 

[u~ (x) , u + (x')] = —iG~(x, x')id, 

[u(x),u(x')} = —iG(x,x')id, 

where id is the identity operator, and the pointwise commutators are interpreted as formally 
expressing the corresponding distributional identity, e.g. 

Hf),u(h)] = -i f f(x)G(x,x')h(x')^^\{x)^-\g\(x') dxdx', 
Jnxn 

where / and g are test functions. 

III. A theory of creation and annihilation operators can be constructed from these formulas 
and (3.7). All the usual propagators can be deduced from the elementary kernels and Q. 

IV. In Minkovski space, the above decomposition agrees with the usual split of u into negative 
and positive energy solutions via Fourier transform. 

One feature that stands out is the existence of positive energy solutions. This seems to be in 
contradiction with the established fact that in a general spacetime there is no natural notion of 
positive frequency (which is itself a consequence of the lack of uniqueness for the vacuum state 
in such situations). Recall, however, that such a notion is well defined for spacetimes that are 
(i) globally hyperbolic and (ii) stationary. The former property has been implicitly employed 
by addressing the problem solely from a local perspective and considering x, x' G D(S) nil in 
Assumption 3.1. These are not essential: when considering a global point of view, it is rather 
natural to restrict oneself to globally hyperbolic spacetimes, and our lemmas and propositions 
can be extended to this setting (recall that we focused on local constructions because the 
approximation arguments we used are local in nature). 

Although not implicitly assumed, property (ii) was "almost" assumed, in the following 
sense. We have supposed the existence of the kernel Q, but sufficient conditions that imply 
the validity of Assumption 3.1 have not been given. It can be shown, however, that Q in fact 
exists for globally hyperbolic stationary spacetimes (Euclidean at infinity for non-compact 
Cauchy surfaces) [M] (see also [C]). But, in the special case of globally hyperbolic stationary 
spacetimes, the decomposition into positive and negative energy solutions does in fact hold 
true [W]. Since it is not know, however, whether global hyperbolicity and stationarity are also 
necessary conditions for the existence of Q [L3], here we preferred to take the slightly more 
general point of view of postulating the existence of Q itself. 
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